Topological properties of superconducting junctions 
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Motivated by recent developments in the field of one-dimensional topological superconductors, 
we investigate the topological properties of s-matrix of generic superconducting junctions where 
dimension should not play any role. We argue that for a finite junction the s-matrix is always 
topologically trivial. We resolve an apparent contradiction with the previous results by taking into 
account the low-energy resonant poles of s-matrix. Thus no common topological transition occur in 
a finite junction. We reveal a transition of a different kind that concerns the configuration of the 
resonant poles. 
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Superconducting junctions, including 

superconducting-normal (SN) ones where dissipa- 
tive conduction can take place and superconducting- 
superconducting (SS) ones where a discrete spectrum of 
bound Andreev states is formed, have been in focus of 
condensed-matter research for almost fifty years [H 
An indispensable compact approach to superconducting 
junctions employs a scattering matrix that relates 
incoming and outgoing wave amplitudes that obey 
the Bogolyubov-deGennes (BdG) equation 0, Q). The 
beauty and power of this approach stems from its ability 
to incorporate numerous microscopic details in a com- 
pact form of the scattering amplitudes. Straightforward 
extensions permit to include magnetism, spin-orbit 
interaction, non-trivial superconducting pairing 
The s-matrix approach can be easily combined with 
semiclassical treatment of electron transport in the 
framework of a quantum circuit theory Q. 

Recent developments in the field of superconductivity 
require revision of the common assumptions concerning 
the structure and properties of the scattering matrix of a 
superconducting junction. Kitaev in 2000 has suggested 
a model Id p-wave superconductor Q that exhibits a 
topological order. It has been shown recently that the 
same topological order can be realized in more realistic 
systems that combine spin magnetic field Q with strong 
spin-orbit interaction [g, |9( . Similar situation would oc- 
cur in a superconductor on the top of topological insu- 
lator or half-metal [10|. The relevance of these develop- 
ments for generic superconducting junctions is not imme- 
diately relevant. Indeed, the general properties of those 
are not supposed to depend on dimension [11], while 
topological ordering considered is specific for one dimen- 
sion [12j thus suggesting that the topological properties 
are not at all manifested in junctions. However, a num- 
ber of spectacular predictions and device schemes that 
relate the topology and junction properties has appeared 
in the last years. Those include: prediction of so-called 
47r periodic Josephson effect @, [lj, E3] , formulation of a 
criterion for topological transition in terms of reflection 
matrix of a junction [l5| , proposals of topological qubits 



based on majorana bound states [a, Il4ll as well as their 
readout with qubits of different type [l6| . 

Thus motivated, we have performed a topological anal- 
ysis of a general BdG scattering matrix concentrating on 
energy dependence of its eigenvalues. This rather elemen- 
tary analysis shows that i. there are topologically non- 
trivial (TNT) s-matrices characterized by real eigenval- 
ues at zero energy, ii. there are topologically non-trivial 
trajectories (TNTT) in the space of topologically trivial 
(TT) s-matrices, that pass a matrix with real eigenvalues 
at E = odd number of times. 

TNT would correspond to a "topological" SN junction 



171 ] , while TNTT would explain 4-7r-periodicity of Joseph- 
son effect in SS junctions [a El. Albeit the same topolog- 
ical reasoning implies topological triviality of all physical 
s-matrices: there are no TNT neither TNTT. This forms 
a paradox that is resolved by recognizing a potentially 
sharp energy dependence of a s-matrix near zero energy. 
Such energy dependence is due to resonant poles [l8| that 
manifest formation and coupling of zero-energy quasilo- 
calized states. With this, we reconcile the predictions of 
d, [lSjj], show the absence of a common topological tran- 
sition and reveal topological transitions related to the 
resonant poles. 

We illustrate these results with two minimal setups, 
SN and SS junctions (Fig. [I}, where a single- channel wire 
with strong spin-orbit coupling and subject to magnetic 
field is brought in contact with a bulk superconductor. 
The Hamiltonian description of this situation is found in 
[8| . In distinction from [8| , we assume finite length of the 
contact. The solutions of BdG equation for a single chan- 
nel encompass spin and electron- hole degree of freedom so 
that the minimal single-channel scattering matrix is 4 x 4. 
The parameter space of the model that includes the su- 
perconducting gap, chemical potential, strength of spin- 
orbit interaction, and magnetic field, can be separated 
into two ranges: "topological" and non-" topological" . 

Let us consider a general s-matrix of a SN junction 
assuming no symmetries. The only constraint on such 
matrix stems from the structure of BdG equation: its 
Hamiltonian satisfies H* = —tiHt\, where the operator 
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FIG. 1: Setups to illustrate general topological properties of 
BdG s-matrices. a. Finite-length wire with strong spin-orbit 
coupling on the top of superconducting lead forming a SN 
junction, b. Finite-length wire between two superconductors 
forming SS junction. Grey ellipses indicate "buried" zero- 
energy states. 
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FIG. 2: Energy dependence of s-matrix eigenvalues, (a) 
Topologically non-trivial (TNT) case, corresponding to the 
"topological" parameter range in (b) Generic topologi- 
cally trivial (TT) case. (Numerical results for the setup in 
Fig. la in the limit L — > oo.) 



Ti switches electrons and holes. The constraint is conve- 
nient to represent in so-called Majorana basis [19J where 
the Hamiltonian is antisymmetric and the scattering ma- 
trix satisfies S(E) = S*(—E) , E being energy counted 
from the chemical potential of the superconductor. Wc 
will consider only energies E within the energy gap of the 
bulk superconductor. In this case, there are no scatter- 
ing waves in the bulk of superconductor, the matrix S is 
in the basis of normal-metal scattering waves satisfying 
unitary condition. 

Let us concentrate on (continuous) energy dependence 
of the matrix eigenvalues e lx ^ E \ That can be repre- 
sented as a manifold of curves in \ ~ E plane (Fig. 
The BdG constraint implies that if a point (x, E) be- 
longs to the manifold, the inverted point (— \, —E) be- 
longs to it as well. These two points can belong to either 
the same curve or to two distinct curves. In the first 
case, the curve is topologically distinct: it is forced to 
pass either x = or x = a t zero energy. If two 
such curves pass the same point, they can be deformed 
by continuous change of Hamiltonian parameters into a 
pair of trivial curves. However, a single curve is topo- 
logically stable: the fact it passes the point cannot be 
changed by Hamiltonian variations. We note that the 
dimension of the physical s-matrices can be always cho- 
sen even. With all this, all s-matrices can be separated 
onto two classes. Topologically trivial (TT) matrices 
have no topologically distinct curves while topologically 
non-trivial (TNT) have two topologically distinct curves 
passing respectively \ = and x — i -71 " at -E 1 = 0. Indeed, 
at zero energy s-matrices are real forming 0(2N) group. 
TT matrices belong to SO(2N) subgroup of 0(2N), while 
TNT belong to 0(2N)/SO(2N). The matrices from these 
distinct submanifolds cannot be continuously deformed 
into one another: indeed, at E = det (TT) = 1 while 
det(TTVT) = -1. 

This classifies s-matrices of SN junction. An SS junc- 
tion is characterized by a combination of two s-matrices 
(Fig. [3]). The spectrum of Andreev states of the junc- 
tion as function of superconducting phase difference <fi is 



obtained from the equation [J] 

= dot (l - S*) ; S = sie l * T3/2 s 2 e-^ r3/2 , (1) 

T3 being Nambu matrix distinguishing electrons and 
holes. It is instructive to note that the unitary ma- 
trix S(cf>) satisfies the same BdG constraint as an SN 
s-matrix. Therefore, the above topological classification 
applies to SS junctions as well. 

In this respect it is crucial to note another topolog- 
ical property that concerns continuous one-parameter 
closed manifolds of TT matrices (trajectories). Intu- 
itively, eigenvalues of a generic matrix " repel" each other 
and never come together. This applies to BdG matrices 
expect a special situation: E = and real eigenvalues. 
Owing to this peculiarity, a trajectory in matrix space 
can in principle pass a matrix where two eigenvalues, 
say, +1, are the same. It turns out that the trajectories 
of the kind can be separated onto two topological classes 
that differ by parity of the number of passes. (Fig. 0} To 
see the possibility for odd number of passes, let us take 
a closed trajectory with a single pass and concentrate on 
two eigenvectors corresponding to the eigenvalue +1. In 
this situation, if the parameter cycles over the trajectory, 
a given eigenvector is transformed not to itself but rather 
to its orthogonal counterpart, this guarantees the stabil- 
ity of this topologically non-trivial trajectory (TNTT). 

Let us understand the results of[l, HH, EH in terms of 
the above classification. Without going into details, we 
enunciate that TNT s-matrices are realized in the "topo- 
logical" parameter range. The TNTT give the topolog- 
ical explanation of the An Josephson effect described in 
these articles. The trajectory parameter in this case is 
the superconducting phase difference <fi. 

However, similar topological considerations show that 
no physical s-matrix belongs to TNT class, neither any 
closed trajectory in parameter space is a TNTT. To prove 
this, let us start with a common (finite) junction mani- 
festing no exotic properties. For our examples, this may 
correspond to zero magnetic field and zero spin-orbit in- 
teraction. The s-matrix at this parameter choice as well 
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FIG. 3: Topological classes of trajectories in the space of TT 
s-matrices. A trajectory is topologically non-trivial (TNTT) 
provided it passes the matrix with two degenerate real eigen- 
values odd number of times. Illustration: the dependencies 
of eigenvalues of the scattering matrix characterizing the SS 
junction on superconducting phase difference <f) for (a) non- 
" topological" and (b) "topological" parameter ranges. 
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FIG. 4: (a): Energy dependence of eigenvalues for NS junc- 
tion in a narrow energy interval illustrates the topological 
triviality of s-matrix for finite length of the contact (L=7 
in units of Dashed lines: "high" -energy TNT eigenval- 
ues. We see the reconnection of neighboring eigenvalues, (b): 
Andreev levels in SS junction versus superconducting phase 
difference at L\ — L2 = 7 (solid lines) as compared to TNTT 
case at LI, 2 = 00 (dashed lines), (c) Energy dependence of 
eigenvalues for case (b) and <fi = n. Dashed lines: TNTT case. 



as all trajectories are topologically trivial. Since there 
is no continuous way to tune scattering matrix from 
TT to TNT class, the s-matrix will stay trivial at any 
strength of magnetic field/spin-orbit interaction, even 
in the "topological" parameter range. This proof is in 
apparent contradiction with the predictions mentioned 



[g, llif , this forming a paradox that motivated us for the 
present research. 

Prior to presenting the solution of the paradox, let us 
mention that the absence of TNTT resolves an annoying 
problem that concerns the parity of particle number of 
the ground state of the SS junction. The level crossings 
at E = are known in the context of ferromagnetic SS 
junctions. Upon passing the crossing, it becomes ener- 
getically favourable to put a single polarized quasiparticle 
to the junction [20]. Therefore, the parity of the ground 
state must be different at two sides of the crossing. Odd 
number of crossings at a closed curve implies indefinite 
parity of the ground state: a situation that is annoyingly 
difficult to comprehend. 

To see how the paradox is resolved, let us consider nu- 
merical results for a finite SN junction in "topological" 
parameter range. (Fig. 4a) If the results are plotted at 
energy scale of the superconducting gap, the pattern of 
energy dependent eigenvalues is apparently of TNT type 
as in Fig. 2a. However, replotting the results near £7 = 
at smaller scale reveals topological triviality (cf. Fig. 4a 
and Fig. 2b). The eigenvalues move fast in the vicin- 
ity of E = reconnecting the branches visible at larger 
energy scale in a rather unexpected way. The typical 
energy scale of such reconnection is small depending ex- 
ponentially on the contact length L, and shrinks to zero 
at L — > 00. This solves the paradox. 

The adequate description of the situation combines a 
smooth energy dependence of s-matrix at E ~ A with a 
pole or poles that are anomalously close to E = 0. Let 
us consider a single pole. The BdG constraints restrict 
it to purely imaginary energy, — iT <C A. The s-matrix 
reads 



S = 1 - 



e-iT 



iT 



-1 |*x*| So, 



(2) 



where is the eigenvector associated with the reso- 
nant level and Sq is the matrix, with smooth energy 
dependence to disregard at E ~ T. The eigenval- 
ues in this energy range are determined from equation 



e/T = E k \^k\ 2 cot(xk 



xi 0) ), exp(ixi 0) ) 



being "high- 



energy" (| E I 3> r) eigenvalues of Sq . They follow the 
pattern in Fig. 2] connecting neighboring "high-energy" 

eigenvalues, exp(ixf ) — > exp(ixjK_x)- This guarantees 
that the total shift of phases of all eigenvalues upon 
crossing a single pole equals 27r. Physically, the pole 
is associated with a quasi-localized zero-energy state be- 
ing formed at the far end of the wire. If the contact 
length exceeds the localization length, this state is ef- 
ficiently "buried" (r <C A) in the superconductor and 
hardly accessible for incoming electron or hole waves ex- 
cept E — when the scattering of the waves become res- 
onant. Andreev conductance of the junction is expressed 
as Ga — GqTx (T3ST3S'). In the resonant energy inter- 
val, the energy dependence of the conductance assumes 
a universal form Ga(E) = Ga + E -F +r -2 (Ga{0) — Ga) 7 
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FIG. 5: (a): Fork SN junction to illustrate topological tran- 
sitions concerning the resonant poles. (b,c): Configurations 
of the resonant poles in the complex energy plane (b) before 
and (c) after a transition. At the transition point, the poles 
are degenerate (double gray circle) 

Ga(0),Ga being its values at E = 0, \E\ > V that de- 
pend on details of the junction. 

Let us turn to the SS junction in the "topological" 
parameter range. Solving Eq. [T] gives the spectrum of 
Andreev states (Fig. [5)3). We observe the level crossing 
at E = 0, 4> = ft being lifted in a narrow energy interval. 
Strikingly, we observe another pair of levels with energies 
remaining small in the whole range of phase. These levels 
are absent in TNTT picture and emerge as a consequence 
of topological triviality of the s-matrix. Since there is no 
level crossing at E = 0, the parity of the ground state is 
always even. 

The situation can be comprehended if we notice that 
each matrix §i,§2 forming the resulting s brings a reso- 
nant pole corresponding to a "buried" zero-energy state 
at far end of each wire. The s thus has two resonant 
poles. The mixing of the two "buried" states results in 
their (phase-dependent) energy splitting and formation 
of the pair of low-energy Andreev levels. The eigenval- 
ues of s-matrix move in the narrow energy interval re- 
connecting next-to-nearest (two poles) neighbour "high- 
energy" eigenvalues (FigJ3t>). This brings four rather 
than two states in the vicinity of the crossing point 
E = 0, <f> = 7T, x = 0) all being involved in the lifting 
of the degeneracy. The detailed theory of the crossing 
point will be presented elsewhere. 

Since the s-matrix remains topologically trivial, there 
can be no sharp transition in its characteristics that 
would correspond to the "topological" transition in the 
(rather unphysical) limit of infinite wire. However, a 
BdG s-matrix with resonant poles is characterized by a 
topological number that can change sharply upon chang- 
ing the parameters. 

Let us illustrate this with a two-pole scattering matrix 
correspond to the fork setup in Fig. 5 a. Here the scat- 
tering matrices S± , S2 of fork tines bring a resonant pole 
each. The BdG symmetry leaves two distinct possibili- 
ties for the poles of the total scattering matrix :i. both 
poles lie on the imaginary energy axis {E — —iTi, — il^), 
ii. they form a pair symmetric with respect to reflection 
ReE — > —ReE (E = ±e — iT). One can now change the 
s-matrix Sq describing the normal scattering in the fork. 



If the tines are open to the lead states, the pole config- 
uration should be like one for two parallel SN junctions: 
the possibility ii is realized. If the tines are isolated, the 
"buried" states mix resulting in an energy spitting: the 
possibility i is realized. We thus expect and [2l[ prove 
the transition at intermediate coupling. 

Generally, one can characterize a BdG s-matrix of ar- 
bitrary dimension with a topological number that is just 
the number of poles lying precisely on the imaginary axis. 
We expect this number to change by 2 upon changing 
the parameters, this gives a series of "topological" tran- 
sitions. (Fig. 5 b,c) Two poles are degenerate at the 
transition point. However, since in general the degener- 
ate poles are at finite imaginary energy T, the manifes- 
tations of the transitions in transport properties are lim- 
ited. The energy-dependent Andreev conductance does 
not seem to have a singularity at the transition point. 

We have performed the topological analysis of the 
properties of SN and SS junctions characterized by BdG 
s-matrices. We have proven topological triviality of phys- 
ical matrices that describe finite-size junctions: there is 
neither TNT, nor TNTT. This implies the absence of a 
sharp "topological" transition upon crossing to "topo- 
logical" parameter range as well as the absence of Air- 
periodic Josephson effect. We have resolved the appar- 
ent contradiction with results of [s|,[l3|,[lH by considering 
the low-energy poles of s-matrices. The resulting sharp 
energy dependence at E Olcads to Lorentian energy 
dependence of Andreev conductance. We have demon- 
strated a topological transition (or a series of transitions) 
of a different kind associated with a change of the con- 
figuration of the resonant poles in complex energy plane. 

This research was supported by the Dutch Science 
Foundation NWO/FOM. The authors are indebted to C. 
W. J. Beenakker, C. Kane, R. M. Lutchyn, F. von Oppen 
and L. P. Kouwenhoven for useful discussions. 
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APPENDIX 

We consider the fork setup presented in Fig. 5 a of the 
main text. Let us specify the normal-scattering matrix 
5o that determines the reflection of electrons and holes 
coming from the lead and their transmission to the fork 
tines. In general form, we can write this matrix in blocks 
of reflection and transmission matrices: 




FIG. 6: Pole positions versus parameter £ characterizing the 
model setup. At £ = there is no reflection in the fork, 
while at £ = 1 the reflection is complete so that the tines are 
isolated. We chose F2 = 2Fi. 




So 



Rd T u d 

Tdu Ru 



(3) 



FIG. 7: Andreev conductance of the fork setup versus energy 
for a set of parameter values (from the top curve at E — 
to the bottom one): £ = 0.05,0.15,0.24 (corresponds to the 
transition, thick line) 0.35,0.45,0.55, Go = 2e 2 /h. 



with Rd(R u ) being the reflection matrix to the lead (to 
the tines). We denote the scattering matrix of the two 

tines as S = diag |j§i,^2|- The total scattering matrix 
then reads: 



St, 



Rd + TduSTud + Td u SR u ST u d 



Rd + Td u S 



1 



1 — R U S 



j-Tud 



(4) 



The poles of S to t are determined by the zeros of the de- 
terminant in the above expression, 



det (I - R U S) = 0. 



(5) 



Now we implement the pole decomposition (Eq. 2 of 
the main text) for each Si. Using the fact that R u is 
invertible we obtain: 



det R-\S: 



diag 



1 - 



2iTo 



2iT 1 



>< -01 1 



E + iT- 



#2 >< 4>2 



= 0, (6) 



with S p being the energy-independent (at E ~ T) part 
of the s-matrix of the tines. We use convenient notations 
E = ie and invert the matrix R^Sp 1 — 1 to arrive at 



x diag 



2Ti 



l^i >< "01 



det ! 1 

2r 2 



SpR u 



1 SpR u 



102 >< lfa\ 



= 0.(7) 



The condition of zero determinant can be then expressed 
in terms of the matrix elements in the basis \ip2 > , |0i > 
that we arrange into a 2 x 2 matrix M: 



M = 



Vi 

1p2 



SpRu 
1 — SpR u 



1p2 




M 



'112- 



2ri 



1 + M 22 



e+r 2 



= 0. 



(8) 



(9) 
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This reduces to a quadratic equation 



e 2 + e(Ti + 2T 1 M 11 + 

+r 1 r a ((i + 2Mii)(i + 2M a2 ) 



r 2 + 21W22) + 

- 4M 12 M 21 ) = 0. (10) 



We note that all matrix elements are real owing to the 
BdG symmetry. This is why the equation roots are ei- 
ther real or mutually conjugated. The intermediate sit- 
uation between these two cases is two degenerate real 
roots. This is the point of topological transition. 

Let us prove that the topological transition of this kind 
should inevitably occur upon changing the So from full 
transmission to full reflection. 

In case of full transmission. R u — > 0. Therefore, M —> 
and the roots e — > T\ or — > r 2 . Therefore, we are 
in the situation with two real roots. The case of full 
transmission is a bit more difficult to handle. In this 
case, R u is a unitary matrix. With this, 



1 



= h - cot(h) 

l-S p R u 2 2 v 1 

h being a Hermitian matrix. The BdG symmetry re- 
quires h to be asymmetric in Majorana basis. Since 
the pole eigenvectors |V>i), IV'2) are real in this basis, 
Mn = M 22 = -1/2 and M 12 = -M 2i . With this, the 
quadratic equation reduces to 



0. 



(11) 



that has two conjugated (and purely imaginary) roots. 

This implies that changing from full transmission to 
full reflection requires passing a point when two roots 
are degenerate, so that, the point of the topological tran- 
sition. 

We illustrate with a simple model involving 4 x 4 s- 
matrices. We choose 



Ru — —Rd = 



1 1 



V2 V-i 1 



T, 



d — Tdu 



Ru- 



(12) 
(13) 



Since these matrices describe normal scattering, they are 
diagonal in electron-hole space. The parameter value 
£ = 0(£ = 1) corresponds to the case of full transmis- 
sion (reflection) We choose the scattering matrix of the 
tines to be 



E 



s,_, 





gl 
E+ir! 







E 



E+iT 2 



»r 2 

E+iT 2 



E+iT-t 



E 
E+iV t 







iT 3 



E+iT 2 



E 
E+iT 2 



(14) 



Here, two upper rows are for electron part of the wave 
while two lower rows are for hole part. This is the sim- 
plest matrix with two poles at — iTi,— iT^. In Fig. [S] 
we plot the positions of poles versus £. The topological 
transition takes place at £ « 0.24. 

We plot the energy dependence of Andreev conduc- 
tance for this model setup at several values of £ in Fig. 
[7] Qualitatively, one expects a single-peak energy de- 
pendence in the case of big transmission, and a double- 
peak dependence in the case of low transmission, the po- 
sitions of the peaks corresponding to the energy of An- 
dreev bound state. This is indeed seen in the Figure. A 
naive expectation would be that the intermediate situa- 
tion where the second derivative of Andreev conductance 
vanishes, occurs at the topological transition. Yet this 
does not happen. 

It remains unclear at the moment whether the topo- 
logical transition under consideration manifests itself as a 
singularity of any physical quantity. While this is likely 
the case in the model of non-interacting electrons, the 
interactions may change this. 



